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A model for cooperative sequential adsorption that incorporates nearest-neighbor exclusion and 
next-nearest neighbor interaction is presented. It is analyzed for the case of one-dimensional dimer 
and two-dimensional monomer adsorption. Analytic solutions found for certain values of the in- 
teraction strength are used to investigate jamming coverage and temporal approach to jamming in 
the one-dimensional case. In two dimensions, the series expansion of the coverage 6{t) is presented 
and employed to provide estimates for the jamming coverage as a function of interaction strength. 
These estimates are supported by Monte Carlo simulation results. 
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Surface adsorption phenomena are important in a 
great number of physical, chemical and biological sys- 
tems. Equally large is the number of phenomena them- 
selves that occur when particles/molecules are adsorbed 
to a surface. The surface-adsorbate interactions can be 
broadly classified into two categories, physisorption and 
chemisorption [1]. While the former is associated with 
electrostatic interactions, van der Waals forces @, Q, 
the latter is commonly associated with the formation of 
chemical bonds. Due to the strength of the chemical 
bond, chemisorption is often irreversible for temperatures 
of interest. An well-known example is the adsorption of 
water on Fe(001) [3]. A simple but effective model for 
such irreversible adsorption is the random or cooperative 
sequential adsorption (CSA) Q: particles are adsorbed 
randomly in a sequential manner without diffusing or des- 
orbing. 

For irreversible adsorption, a central quantity of inter- 
est is the coverage, 9, of the final adsorbed monolayer 
attains, the jamming coverage, 8j. Detailed knowledge 
of this jamming coverage might, for example, become 
important for chemical sensing devices, such as micro- 
cantilevers Q, in order to distinguish between different 
species of adsorbates. The jamming coverages have been 
estimated within several approaches [5J. However, impor- 
tant and significant effects due to interactions between 
adsorbates on the jamming coverage have not been con- 
sidered in detail and are the focus of this work. 

One of the simplest models for interaction between ad- 
sorbates is the nearest-neighbor exclusion (NE) mecha- 
nism that causes an adsorbed particle to block its nearest- 
neighbor binding sites from adsorption (4j. It has been 
shown 7] that for the NE process on a two-dimensional 
surface with a square-lattice arrangement of binding sites 
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the jamming coverage approaches a non-trivial value of 
9 j = 0.3641 which is lower than that of the ideal limit, 
#o = 0.5. This is a consequence of the stochastic nature 
of the adsorption process which results in the exponen- 
tially small probability of an ideally covered surface. 

Here, we introduce, in addition to the NE, a short- 
range adsorbate-adsorbate interaction which provides a 
more accurate and realistic description of adsorption pro- 
cess. This interaction affects the rate of adsorption of a 
next-nearest neighbor (NNN) binding site and influences 
the jamming coverage. Such interactions might be caused 
by a variety of mechanisms. If the adsorbates have an ef- 
fective charge, then resulting electrostatic forces would 
lead to a repulsive interaction: for example, hydrogen on 
Pd(100) acquires a dipole moment due to charge trans- 
fer from the surface In this case, assuming that the 
adsorption rates are of Arrhenius type and that each oc- 
cupied NNN of an available binding site contributes an 
equal amount to the binding energy e- m t, the adsorption 
rate including the interaction with n occupied NNN could 
be modeled by r n = ro exp (— ne; nt /T) Q (where ro is 
a typical rate of adsorption and T is the temperature). 
Attractive interaction might arise if an adsorbate (e.g. 
water on Pd(10Q) 0. Q5) 

induces a local change in the 
surface structure [IlJ that increases the rate of adsorp- 
tion at NNN sites. Another mechanism that would lead 
to an effective attractive interaction could involve pre- 
cursor layer diffusion pjj]: a gas particle might become 
physisorbed even if it collides with an already adsorbed 
particle [l^]. In that case, it might either desorb or dif- 
fuse to the next available site surrounding the adsorbed 
cluster for chemisorption. 

Below, we use a linear approximation in n for the ad- 
sorption rate, r n = 1 — ne (ro = 1 by rescaling the time). 
Such a linear dependence on n arises naturally for the at- 
tractive interaction of adsorption via the precursor state, 
as each occupied NNN should contribute equally to the 
adsorption of the surrounded site [HI]. It is clearly the 
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first order approximation of the Arrhenius-type rate pre- 
sented above, which is valid for high temperatures or 
small interaction energies £j n t and also a possible as- 
sumption for the mechanism of morphology changes in 
the surface. It should be mentioned, that the effects of 
adsorption rates that depend on NNN occupation, espe- 
cially on the island structure, have been previously con- 
sidered [lj|, albeit with a different choice of rates. 

The effects of the NNN interaction on adsorption can 
be investigated by means of rate equations Q which 
describe the evolution of the marginal probability den- 
sity P(G; t) of finding a configuration G of lattice sites 
empty at time t, irrespective of the state (occupied or 
vacant) of the remaining sites of the lattice. For the two- 
dimensional square lattice, the rate equations are [Til ] 



d t P(G; «) = - X) S r nP({G U A}„; t) 



(1) 



The only way that in an irreversible adsorption process 
the probability of finding a set G of sites can change is 
by adsorption at one of its binding sites i <G G. Due 
to NE, i must have z empty nearest neighbors. There- 
fore if i lies on the boundary of G, this can only hap- 
pen if G is a subset of the larger set of empty sites, i.e. 
G U Di. The subscript n in {G U -D;}„ refers to the ad- 
ditional interaction with the environment of site i sur- 
rounded by n occupied NNN. Using the fact that the 
marginal probability densities obey the following rela- 
tion, P(G U {aj = 1}; t) + P(G U {cr 3 = 0}; t) = P(G; t) 
(<7i denotes occupation of site i) we can recast the RHS 
of Eq. (fTJ) completely in terms of probability densities of 
configurations of empty sites. For example, considering 
only the contribution to the rate equation due to adsorp- 
tion at the dotted site, Eq. (TTJ) reads 



d t P(o§e;t) 



(l-2e)P(°§§°;i) 



-2eP(o88o;f) 



(2) 



Formally, we write such a rate equation as dtP(G; t) = 
—CP(G;t) where C is the operator that generates the 
configurations G' of empty sites that can produce G by 
a single adsorption event. 

First, we analyze the situation in Id. Monomer ad- 
sorption with NE and NNN interaction is equivalent to 
dimer adsorption with with nearest-neighbour interac- 
tion in Id [Tt], [HI, which we will consider here. This 
model has been solved for general cooperative rates [lj| . 
However, the temporal approach to jamming, which is 
mainly of our interest, is not readily available from such 
a solution. Therefore, we present here a different form of 
the solution that is suitable for our purposes. 

For this process, the following rate equations for find- 
ing a stretch of m vacant sites can be written: 

d t P(l;t) = -2(l-2e)P(2;i)-4eP(3;t) (3) 
d t P(m;t) = — [m — I — 2e]P(m; t) — 2P(m + l;t) 

-2eP(m + 2;t) for m > 2 . (4) 



With the initial condition of an empty lattice, Eq. ([4]) is 
solved exactly by 



P(m;t) = exp[-(m - I - 2e)t - 2(1 - e 



-e(l 
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(5) 



for m > 2. The solution for the case m = 1 (see Eq. 
is given by 

P(l;i) = 1 - 2e- 2 " e ((1 - 2e)P 2 (t) + 2e/|(t)) (6) 
where 



P(m; t') dt' . 



(7) 



The probability P(l; t) is particularly important for eval- 
uation of the critical coverage, 6j = 1 — lim^oo P(l;i). 
The integrals I^it) can be evaluated analytically only 
for some special cases. Namely, by considering negative 
values for e = —a (a > 0), we can write I^if) in terms 
of a sum of lower incomplete gamma functions 7 20] , 
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k=0 



m + 2a~2 
k 



(8) 



where = e 1 /" a 3 ^- m - 2a , w(a,t) = ^(e'* - 1/a). 
The infinite series only converges in the interval a S 
[0,2), but it becomes finite for half-integer and integer 
values a = n/2 with n — 0, 1, 2, . . .. Using the identi- 
ties 7(0 + 1, x) — 07(0, x) — x a e~ x , 7(1, x) = 1 — e~ x and 
7(1/2, x) = y / 7rerf(-y/S) [2(|, we find, for example, P(l; t) 
for e = -1/2 to be 



P(l;t) = 1 + 2^ 
and for e = — 1 to be 
P(l;i) = l- 



^erf(u)) -(2 + V2w), 



u)(l/2,0) 



(9) 



'- erf(io) - (2 + 3 w + 2 w 2 ) e 



w(l/2,t) 



uj(1,0) 

(10) 

w(l,t) 



where erf(x) denotes the error function. It follows from 
Eqs. ([9]) and fTO]) that the time-dependent coverage, 
9{t) = 1 — P(l;t), asymptotically approaches the crit- 
ical value 8j ~ 0.876681 as 



0(t) 



(eri(-V2) - erf(-l/V2)) 



4e -f- 2 * 



0((«;(l/2,i)-i£7 



and 8 



j ~ 0.885296 as 

-1 



6{t) = e- 1 -2 + erf(-l)/2 + 2e" 1 " 3t 
+0((w(l,t)-w oo (l)) 4 ) 



(11) 



(12) 
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for e = — 1/2 and e = — 1, respectively. These expressions 
have been obtained by expanding Eqs. © and (fTU|) about 
Woo (a) = hmt^oo ■u;(a, t) = —l/y/a such that w — = 
^fae~ l . The time-dependent coverage for any e = — n/2 
can be obtained analytically in a similar manner. The 
important point to note here is that Eqs. (fTTj) and (fT2"l) 
suggest that for a given n, all terms of order smaller or 
equal n will drop out of the expansion so that the leading 
time-dependent term in the approach to jamming when 
t — ► oo is exp(— t/r) with t = l/(n+l). Extrapolating to 
any value for e, the characteristic time r to jamming for 
the one-dimensional process is then given by t — 1/(1 — 
2e) which can easily be verified numerically. This form of 
the temporal approach to jamming also follows from the 
time dependence of the sticking probability P(2; t) given 
by Eq. ©. 

The 2d case of monomer adsorption with NE and NNN 
interaction cannot be handled analytically and the meth- 
ods of series expansion (SE) and Monte Carlo simulations 
(MC) were employed. The SE is an Taylor expansion of 
any probability density P(G; t) in t using the rate equa- 
tions ©, P(G;t) = £ n=Q i-^£™P(G, i) |, 0. For 
the coverage, the expansion of P(o\t) is of interest. We 
have implemented the algorithm introduced in (2lj with 
the extension of allowing for polynomial values in e in or- 
der to compute the coefficients of the series numerically. 
The series coefficients up to order 13 were computed (see 
Tab. H| using this algorithm. With these coefficients and 

the first terms in the expansion of P( ogo ; t) = P(G\\t), 

we can also calculate the expansion of the sticking prob- 
ability, S(9) = P(Gi,9), as a function of the coverage 
6{t). The first few terms are 

S{9) = 1 - 5 6» + (6 - 4e) 6» 2 + |(1 + 3e - 5e 2 ) 6> 3 

-|(l-26e- 13e 2 + 84e 3 )6» 4 + 0(6» 5 ) (13) 

The interaction hardly affects the sticking probability for 
small coverages (see inset in Fig. [T]) and the effect can be 
seen only close to the jamming limit. This is what one 
would expect, given the short-range nature of the interac- 
tion. The value of the sticking probability for attractive 
(repulsive) interaction are expectedly greater (smaller) 
than for the adsorption without NNN interaction. 

To obtain estimates for 9j from the series the stan- 
dard Pade approximants @, H3] [n, m] (where n and m 
are the orders of the polynomials in the numerator and 
denominator of the Pade approximant) were used. First, 
the transformation of variables, y = (1 — e xpf — (1 — 
be)t))/(l — be) (similar to that used in Ref. 23]), was 
carried out with b being an adjustable parameter. It 
is clear from the preceeding discussion that this mim- 
ics the approach to jamming in Id where b take the 
value 6=2. It has been found before Q that using 
the knowledge of the temporal behavior of the coverage 
in Id for the transformation of variables can consider- 
ably improve estimates in 2d. However, instead of choos- 



TABLE I: The series coefficients c nm for the Taylor expansion 
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ing a value for 6, we use it to make the estimates of 0j 
independent of the choice of Pade approximant. In or- 
der to do so, the free parameter b has been found by 
minimizing the cost function C(e, b) = [6j[e, [6, 6],&) — 

6j(e,[6,7],b)) 2 + (6j(e,[6,6],b) - 6j(e, [7, 6], b)f + 

(0j(e, [6, 7], b) — 9j(e, [7, 6], b)) with respect to b where 
(?j(e, [n, m], b) is the jamming coverages obtained for the 
highest-order Pade approximants available. 

The results of this analysis are presented in Fig. Q] and 
compared with the values of jamming coverage calcu- 
lated numerically by MC simulations. In the MC sim- 
ulations, an event-driven algorithm 0, [l6|, HH was used. 
Within this algorithm, all the susceptible binding sites 
were grouped depending on the number of occupied NNN 
sites. A binding site for the next adsorption event is then 
drawn randomly out of a group according to the rates r n 
and the waiting times are distributed exponentially. The 
results for the jamming coverage calculated numerically 
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for a wide range of interactions are presented in Figs. [T] 
and [2] 



In Fig. [I] one can clearly see that the series expansion 
and the MC simulations agree for e £ [—2.0,0.1]. The 
relative effect of the interaction is strongest around the 
point e = and then flattens as e becomes more negative, 
i.e. the interaction becomes more attractive. In Fig. [2j 
we can see that only for very large negative e, e < — 10 4 , 
the jamming coverage comes within a few percent of the 
ideal coverage of 9j = 0.5. 



FIG. 1: (color online). Comparison of jamming coverage 8j 
for e = —5.0, . . . , 0.24 from the series expansion up to order 
13 using Pade approximant [6, 6], [6, 7] and [7, 6] for the opti- 
mization described in the text (+) and from MC simulations 
for a 200 x 200 - lattice and 100 iterations (red — ). Errors 
in simulation data points are of order 0.001. Inset: Sticking 
probability S as function of coverage 6 up to order N = 9 
for three different e = 0.1( — ), e = 0.0(— ), and e = -0.5 
(•■--)■ 
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In conclusion, we have presented the analysis of a co- 
operative sequential adsorption model, that takes into 
account the effects of nearest-neighbor exclusion as well 
as physically important (repulsive or attractive) inter- 
actions between next-nearest neighbours. A one- and 
two-dimensional process, dimer adsorption and monomer 
adsorption with nearest-neighbor exclusion, respectively, 
have been studied. For the one-dimensional process, we 
computed the coverage analytically for for a family of spe- 
cial cases, where the interaction parameter e takes nega- 
tive half- and integer values. This allowed us to compute 
the jamming coverage and to extract the temporal ap- 
proach to jamming. For the two-dimensional process, we 
have computed the series expansion for the coverage as 
a function of time and have found the jamming coverage 
for various strengths of interaction. Monte Carlo simu- 
lations convincingly support the series expansion results 
and provide estimates for the jamming coverage that are 
unaccessible to the series expansion, thus demonstrating 
a slow convergence to the ideal coverage with increasing 
attractive interaction. 



FIG. 2: (color online). Jamming coverage 8j for e = 
-1.0, . . . , -1.0 x 10 6 from MC simulations for a 200 x 200 
- lattice and 100 iterations. Errors in simulation data points 
are of order 0.001 and thus smaller than the symbol size. 
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